Large Quantum Superpositions and Interference of Massive Nanometer-Sized Objects 
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We propose a method to prepare and verify spatial quantum superpositions of a nanometer- 
sized object separated by distances of the order of its size. This method provides unprecedented 
bounds for objective collapse models of the wave function by merging techniques and insights from 
cavity quantum optomechanics and matter wave interferometry. An analysis and simulation of 
the experiment is performed taking into account standard sources of decoherence. We provide an 
operational parameter regime using present day and planned technology. 
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Quantum superpositions of a massive object at two 
spatial locations are allowed by quantum mechanics. 
This puzzling prediction has been observed in seminal 
matter wave interferometry experiments with electrons, 
neutrons, atoms and dimers, van der Waals clusters, and 
even complex molecules (e.g. C70, C60F48) [1]- Preparing 
quantum superpositions of even larger objects is consid- 
ered to be extremely challenging due to the decoherence 
caused by interaction with the environment [2]. How- 
ever, succeeding in this task would allow completely new 
tests of quantum mechanics: this includes experiments 
in a hitherto unachieved parameter regime where col- 
lapse theories predict quantum mechanics to fail [3, 4], 
or even more general tests of quantum theory against 
full classes of macrorealistic theories [5]. Moreover, these 
states would be so fragile to environmental interactions 
that one could exploit this ultra-high sensitivity to design 
a new generation of sensors. Pushing large objects to the 
quantum regime is also the aim of cavity quantum op- 
tomechanics [6]. Similarly to laser cooling of atoms, the 
radiation pressure of light is exploited to cool and coher- 
ently manipulate the mechanical motion of some degree 
of freedom (e.g. the center of mass) of a massive object 
and even to create quantum superpositions of harmonic 
vibrational states [7, 8]. 

In this Letter, we present a method to prepare spatial 
quantum superpositions of massive objects (with masses 
of ~ 10 7 amu) based on cavity quantum optomechanics 
and show how it can be used to test wavefunction col- 
lapse models. This builds upon the recent proposal of 
using an optically levitating nano-dielectric as a cavity 
quantum optomechanical system [8-11]. The main idea 
is to trap a dielectric sphere in the standing wave field of 
an optical cavity. The mechanical motion of the sphere's 
center of mass along the cavity axis is predicted to be 
a high-quality mechanical oscillator due to the absence 
of clamping losses. This facilitates laser cooling to its 
motional ground state (see also experiments on feedback 
cooling of an optically levitated microsphere [12]). In ad- 
dition, a cooled levitating object offers the possibility to 
be released by switching off the trap [10], creating in this 
way a scenario similar to matter wave interferometry ex- 
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FIG. 1: (Color online) Schematic representation of the pro- 
posal, (a) The optically trapped object is laser cooled using a 
high- finesse optical cavity, (b) The trap is switched off and the 
wave function expands during some time t\. (c) The object 
enters into a second small cavity where a pulsed (of time r) 
interaction is performed using the quadratic optomechanical 
coupling. The homodyne measurement of the output phase 
measures x 2 and prepares a quantum superposition state con- 
ditional of the outcome result pl- (d) The particle falls during 
a time £2 until its center of mass position is measured, which 
after repetition unveils an interference pattern for each pl- 



periments. Here, we will use precisely this feature both 
to coherently expand the wave function over a large spa- 
tial region and to enhance the non-linear coupling that 
is required to prepare large quantum superpositions. 

More specifically, the linear and quadratic coupling in 
cavity optomechanics after displacing the cavity field (see 
e.g. Sec. V.A.I and App. B.2 in [10]) is given by 

Hom = —fig (a + cJ)x + hg q (a + a))x 2 , (1) 

where a(a)) is the annihilation (creation) operator of 
a cavity photon, x = x/xq is the dimensionless posi- 
tion operator of the mechanical resonator, with xo — 
[h/(2muj t )} 1 ^ 2 its zero point motion, m the mass, and uj t 
the mechanical frequency. The photon-enhanced linear 
optomechanical coupling is given by g and the typical 
quadratic coupling by g q = k c x§g^ where k c = 2tt/X c is 
the wave number of the cavity mode. When the equilib- 
rium position of the mechanical oscillator is at the node 
(anti node) of the standing wave, g 7^ and g q = 
(g = and g q 7^ 0). A fundamental figure of merit 
of the cavity-mechanical system is the cooperativity pa- 
rameter defined as C\ = g 2 /(nY) for the linear coupling, 
and C q = g 2 /(nY) = C\ x (k c xo) 2 for the quadratic one. 
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Here, k is the decay rate of the cavity field and T the de- 
coherence rate of the mechanical motion. Ground-state 
cooling requires C\ > 1, whereas non-linear effects, such 
as energy quantization detection [13] or preparation of 
non- Gaussian states without using hybrid systems or sin- 
gle photon resources, require C q > 1. The latter is a very 
demanding condition due to the strong reduction given 
by (k c xo) 2 <C 1. In this Letter we propose to achieve this 
challenging regime by expanding the wave function to a 
given size (x 2 ) ~ a 2 ^> xq, such that 



C q = ^=dx(k c a) 2 , 



(2) 



where g q and f are defined below. Thus, for sufficiently 
large a and Cj, the non-linear regime C q > 1 can be at- 
tained. We remark that this technique is also applicable 
to other setups where the mechanical frequency can be 
varied and hence the wave function of the mechanical 
oscillator expanded [14, 15]. 

We discuss now the different stages of the proto- 
col using levitated nano-spheres [8-11] trapped within 
an optical cavity (Fig. la). The optomechanical cou- 



pling is given by g 



/n~^e c klcV/{W c ), where 



3Re [(e r — l)/(e r + 2)] depends on the relative dielec- 
tric constant e r , c is the vacuum speed of light, V 
the volume for a sphere of radius R and mass m, V c 
the cavity volume, and n p h the cavity photon num- 
ber in the steady state. The decoherence rate of the 
center of mass motion is dominated by light scattering 



with a localization rate 
10]. The decay rate of 



and is given by T = A sc Xq, 
A sc = eln ph cV 2 k 6 J(QirV c ) [9, 
the cavity also has a contribution due to light scatter- 
ing given by ft sc = e 2 V 2 k^c/ (16ttV c ). Sideband cooling 
of the center of mass motion allows the preparation of 
thermal states to a final number occupation given by 
n ~ [k, I (4:Uj t )] 2 +C Z _1 [16], where backaction heating con- 
tributes with Ci ~ c/ (k 2 V c n) for a levitated object at very 
low pressures. Moderate cooling along the other direc- 
tions is also applied to keep the trap stable at low pres- 
sures and reduce the position fluctuations during the time 
of flight. After cooling, the harmonic trap is switched off, 
the object falls (see Fig. lb), and the state evolves freely 
according to 
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The position-localization dissipation part of this master 
equation describes standard decoherence processes such 
as scattering of air molecules and emission, absorption, 
and scattering of black body radiation [2] with the total 
localization rate A s d given below. We remark that de- 
coherence due to light scattering is absent during the 
time of flight since the lasers are switched off. Since 
both the initial state and the master equation Eq. (3) are 
Gaussian, the evolved density matrix can be fully deter- 
mined by computing the moments (x 2 (t)), (p 2 (t)), and 



(x(t)p(t)), where (x 2 (0)} = (2n + l)x% and (p 2 (0)) = 
(2n + 1)^ 2 /(4xq). The spatial coherence length £/, ob- 
tained by noticing that (— x/2\p\x/2) oc exp [— x 2 /^ 2 ], is 
given by £ 2 = 8(x 2 )(V) 2 , where (V) is the mean value of 
the parity operator. 

After an expansion of duration ti, a second cav- 
ity is used to implement an optomechanical double slit 
(Fig. lc). To this end, the setup is aligned such that the 
object passes through a small high-finesse optical cavity 
at an antinode of the cavity mode. Simultaneously, a 
pulse of length r w 2i\ j k is fed into the cavity such that 
a short interaction is triggered. Note that during this in- 
teraction, standard decoherence and, in particular, light 
scattering decoheres the state of the system with a rate 
given by f = A sc (x 2 (ti)) . This can be taken into account 
by adding the corresponding contribution of time r to the 
moments of the Gaussian state before the measurement. 

Linear pulsed optomechanics has been discussed re- 
cently for tomography and cooling applications [17]. 
Here, we extend these results to the case of the quadratic 
coupling (see also [18]). The interaction Hamiltonian, in 
the displaced frame and in the rotating frame with the 
resonant laser frequency, is given by H = p 2 /(2m) + 
hgqy/riphX 2 + hg q (a) + a)x 2 . A key remark is that, at 
this stage, the dimensionless position operator is defined 
as x = x/cr(ti) (hereafter we define a 2 = <J 2 (ti) = 
Xq + h 2 t 2 /(4xQm 2 )). Then g q = g q (a/xo) 2 is the 
quadratic optomechanical coupling enhanced by the en- 
larged wavefunction. Note that, as mentioned above, 
the kinetic term can be neglected since r(p 2 ) / (2mh) w 
uo t T/A <C 1 for short cavities where n ^> uo t . The 
squared position measurement is performed by measur- 
ing the integrated output quadrature of the light field 
Pl = J r ^[^out(^) + ^out(^)]/v^ r - Using the input-output 
formalism, a out (t) -\- di n (t) = \/2^(a + y^h), one obtains 
that (p L ) = X (x 2 ) and (f L -(p L ) 2 ) = l/2+ X 2 (x 4 -(x 2 ) 2 ) 
(we assume a coherent drive such that the optical in- 
put phase noise is 1/2). Therefore, the measurement 
strength of the squared position measurement is defined 
as x — '2\fCq (the physical parameters are chosen such 
that r w 1/r « 27r/ft, see below). Note that the mea- 
surement strength is intimately related to the enhanced 
non-linear cooperativity, see Eq. (2). The generalized 
measurement operator for the measurement outcome pl 
of the integrated optical phase pL is given by 



M = exp 
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where <\> = g q ^/n^r is the phase accumulated during 
the interaction with the classical field. The density ma- 
trix after the measurement is described by pfa + r) = 
AipM. ^ /tr[A4pAi^]. The action of the measurement op- 
erator, Eq. (4), is to prepare a superposition of two 
wave packets separated by a distance d = 2cj^/pl/Xi 
and a width given by approximately a 2 ~ cr/ {^Plx) = 
a 2 /(2dx)- This can be intuitively understood as a con- 
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sequence of the projective nature of the pulsed measure- 
ment [17]: for the ideal case, this measurement prepares 
the system in an eigenstate of the x 2 operator, which 
for a pure initial state with even parity is of the type 
\x) + | — x), i.e. a coherent spatial superposition. The 
separation of the wave packets, determined by the 
outcome of the measurement, represents the effective 
slit separation. In order to prepare (with high proba- 
bility) and resolve the peaks of the superposition state, 
one requires a > d > 2&2 • This sets up an upper 
bound d^ax = a an d a lower bound d m i n = v^/2/x 
for d. A second upper bound is provided by the deco- 
herence during the expansion of the wave function; we 
demand d < dj^ ax = £/• Finally, the total number of 
photons n p h used in the pulse and the time of flight t\ 
are fixed by enforcing that <p compensates the complex 
phase accumulated during the time of flight, which is 
given by ~ (x(t\)p{t\) -\- p(ti)x(ti)) / (4h), as well as by 
fulfilling the condition r w 1/f « 2tt / 'k. This corre- 
sponds to choosing n p h ~ (2n + 1) /[327rCi(k c xo) 2 ] and 
t\ w 16^C l k 2 /[uj 2 (2n + l) 2 A sc /n ph ]. 

After the preparation of the superposition state by the 
pulsed interaction, the particle falls freely during another 
time of flight of duration t^. An interference pattern in 
the mean value of the position is formed with fringes sep- 
arated by a distance Xf = 27rht2/(md). The final step 
of the protocol is thus to perform a position measure- 
ment of the center of mass (Fig. 2d). This requires a 
resolution Sx < Xf, providing a third upper bound for 
d, d^ax — 27rht2/(m5x). Note that sufficiently long time 
ti ~ ma 2 /(hx) is needed to guarantee the overlap of the 
two wave packets. The effect of standard decoherence on 
the visibility of the interference pattern can be obtained 
by solving the evolution of the position distribution for a 
non-Gaussian state under the evolution of Eq.(3). This 
is given by the closed expression [20] 

r°° e -y 2 /° 2 b (t) 
(x\p(t)\x) = / dy — — (x + y\p A=0 (t)\x + y), (5) 

where pA=o(t) is the state obtained with the evolu- 
tion due to the Schrodinger equation only, that is with 
A = 0. As observed in (5), the effect of decoherence 
is to blur the position distributio n with a blurring coef- 
ficient given by a b (t) = 2hm- 1 ^tlA/3. Therefore, the 
fringes separated by a distance x f will be visible provided 
Xf > cr&(£ 2 )/2, which provides the fourth upper bound 
^max = 7r\/ 3 /( t 2Asd)/2. Putting everything together, 
the operational regime for the experiment modelled here 
is given by d min < d < min {< ax , < ax , d^ ax , c^ ax }. 

We now address the experimental conditions re- 
quired for this experiment. The localization rate for 
black body radiation Abb has contributions due to 
scattering A b b, sc oc R 6 T®Re [(e hh - l)/(e b b + 2)] 2 , and 
emission(absorption) of blackbody radiation A b b, e (a) oc 
R 3 T^ e) lm [(e b b - l)/(ebb + 2)], see [2, 9] for the exact ex- 
pressions, ebb is the average relative permittivity, which 
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FIG. 2: (Color online) The operational parameter regime for 
the optomechanical double slit distance d and the the diam- 
eter of the sphere D is plotted (see legend for the lower and 
upper bounds). The simulation of the interference pattern 
is computed for a sphere of 40 nm and d = 0.7 D (circle), 
d = D (triangle), and d = 1.3D (square) in units of D. 
The solid blue (dashed grey) line is the simulated interfer- 
ence pattern with (without) standard decoherence. The dot- 
ted red line is the interference pattern in the presence of the 
CSL model with A = 10 4 Ao s _1 (the upper bound <imax m 
the operational parameter plot provided by the CSL model is 
also shown, see legend). Experimental parameters for the 
environmental conditions: P = 10 -16 Torr, T e = 4.5 K, 
Im[(e b b - l)(e b b + 2)] = 0.1, Re[e bb ] = 2.3, ft = 0.1; for the 
cavity: finesse T — 1.3 x 10 5 , length 2 /xm, waist = 1.5 //m, 
A c = 1064 nm; and for a silica sphere: e r = 2.1 + i2.5 x 10 -10 , 
density = 2201 Kg/m 3 , u t /27T = 135 KHz, and Sx = 10 nm. 
Using this, for a sphere of 40 nm and slit length d — D, one 
obtains k/2tt = r" 1 = 2.8 x 10 8 Hz, Ci = 1500, n ph = 272, 
Ti 206 K, ti 3.3 ms, t 2 125 ms, and a/x 2928. 

is assumed to be constant across the relevant blackbody 
spectrum, and T^ e ) is the internal (environmental) tem- 
perature. Ti at very low pressure can be computed using 
the balance between the emitted blackbody power and 
the light absorption during the optical cooling and trap- 
ping [9]. Second, decoherence due to air molecules is 
described by the master equation Eq. (3) [26], with the 
parameter given by [2] A a i r = 8V2nm a vPR 2 /(3\/3ft 2 ), 
where P is the air pressure, m a is the mass of the air 
molecules and v their thermal velocity. The total stan- 
dard decoherence rate is thus given by A s d = Abb + A a i r . 
The overall performance of this challenging experiment 
is mainly limited by the quality of the cavity used in 
the measurement and the vacuum and temperature con- 
ditions required for the the environment. In particular, 
very good vacuum conditions are needed to keep the co- 
herence of these fragile states. Note however that pres- 
sures down to 10 -17 Torr at cryogenic temperatures of 
T = 4.5 K were reported in [22]. Extremely good cav- 
ities are needed in order to obtain a large cooperativity 
Ci, for instance, consider fiber-based Fabry-Perot cavities 
of length of 2 pm and finesse T ~ 1.3 x 10 5 as discussed 



in [21]. In Fig. 2 the operational parameter regime is 
shown for different sphere sizes and superposition dis- 
tances with the particular set of experimental parameters 
given in the caption. The interference pattern simulated 
by solving the master equation numerically, which de- 
scribes the evolution of the state during the experiment, 
is also plotted. Spheres of ~ 40 nm with a mass of ~ 10 7 
amu can be prepared in a superposition of locations sepa- 
rated by a distance equal to their diameter. In principle, 
the scheme can be applied to even larger objects albeit 
with further constraints on the experimental parameters. 

To conclude, we shall discuss the application of us- 
ing this experiment to test theories beyond quantum me- 
chanics that provide an objective collapse of the wave- 
function for sufficiently large objects. In particular, we 
focus on the paradigmatic model associated to Ghirardi- 
Rimini- Weber-Pear le, see [4] and references therein, de- 
noted as the continuous spontaneous localization model 
(CSL). This theory is derived by adding a non-linear 
stochastic term to the Schrodinger equation. The model 
recovers all the phenomenology of quantum mechanics for 
elementary particles but predicts a fast localization (col- 
lapse) of the wavefunction for larger objects. This comes 
at the price of introducing two phenomenological con- 
stants given by or 1 ! 2 « 10 _T m (related to the localiza- 
tion extension) and Ao ~ 2.2 x 10 _1T s _1 (related to the 
intensity of the localization). For a spherical body [23], 
the CSL model can be cast into a master equation of 
the form of Eq. (3) with A C sl = m 2 \ af( y /aR)/(2m 2 ) ), 
where mo is the mass of a nucleon, and the function f(x) 
defined in [23] has the following limits: /(l) ~ 0.62, 
f(x <C 1) = 1, and f(x ^> 1) « 6x -4 . Recently, 
Adler [24] has reexamined the CSL theory and, by consid- 
ering the collapse of the wave function at the latent image 
formation level, he predicted a significantly larger value 
for Ao, namely Aa = 2 x 10 9±2 Ao- This prediction cannot 
be tested by current experiments [3]. In Fig. 2 we show 
however that a possible CSL process would have a strong 
impact on our experiment already for A = 10 4 Ao (see the 
upper bound c^ax provided by the blurring of the inter- 
ference pattern). The effect is also clearly visible in the 
simulation of the interference pattern. Thus, the experi- 
ment proposed here puts unprecedented bounds for one of 
the most studied collapse models and even challenges the 
recent theoretical prediction given by Adler [27]. Finally 
we note that our scheme allows to prepare superpositions 
of macroscopically distinct spatial states of a massive ob- 
ject. In combination with the specific time-of- flight evo- 
lution this may provide a rigorous experimental test of 
some of the crucial assumptions of macrorealism [5]. 
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